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COMPUTATIONAL  METHODS  FOR  LOW  VELOCITY 


PERFORATION  OF  METALLIC  PLATES 


1 .  INTRODUCTION 


The  ability  of  a  projectile  to  perforate  a  target  can  be  assessed  by 
comparing  Its  kinetic  energy  with  the  mechanical  work  required  to  produce  a 
hole  of  the  appropriate  size  and  shape  in  the  target.  Estimation  of  the 
mechanical  work  to  produce  a  hole  requires  a  knowledge  of  the  deformation 
mechanisms  and  depends  in  a  complex  way  on  such  aspects  as  shape,  size, 
velocity  and  mass  of  the  penetrator,  and  the  impact  obliquity,  as  well  as  the 
density,  mechanical  and  thermodynamic  properties  of  the  target  material.  Two 
basic  approaches  have  been  adopted  towards  the  solution  to  this  type  of 
problem,  viz.  analytical  models  and  computer  code  techniques. 


Analytical  models  of  penetration  are  developed  by  limiting  the  scope 
of  the  problem  to  narrow  geometric  conditions.  An  early  example  was  the  work 
of  Taylor  (1)  on  the  perforation  of  sietalllc  plates  by  pointed  projectiles, 
followed  notably  by  Thomson  121  who  produced  a  solution  for  dishing  failure  of 
thin  plates.  A  number  of  models  of  other  failure  modes  have  been  produced 
(3-51,  which  all  have  simple  analytical  solutions  with  an  equation  relating 
projectile  critical  velocity  to  target  dimensions  and  material  properties, 
in  an  effort  to  better  describe  the  sequence  of  events  during  penetration  more 
complex  analytical  approaches  have  developed  out  of  the  simpler  solutions,  in 
particular  looking  at  plugging  failures  which  are  produced  by  flat  ended 
projectiles  t<,71.  A  number  of  classifications  «f  target  failure  have  been 
produced  and  these  are  presented,  with  discussion  of  concepts  and  models  in  a 
number  of  reviews  I  a— l o 1 .  The  analytical  approach  suffers  from  a  loss  of 
generality,  however  it  has  several  advantages  in  allowing  efficient  parametric 
studies,  demonstrating  the  physics  clearly  and  highlighting  important  aspects, 
and  the  techniques  generally  allow  some  solution  to  a  problem  within 
reasonable  time  and  cost  constraints.  some  analytical  solutions  are 
predictive  because  assumptions  on  hole  size  and  deformation  mode  built  into 
the  model  allow  a  result  to  be  obtained  with  limited  impact  configuration  and 
material  property  data.  other  approaches  are  descriptive  as  they  allow  the 
progress  of  perforation  to  be  analyzed  but  only  after  soote  post  perforation 
measurements  on  hole  geootetry  have  been  made. 
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The  second  approach  is  to  use  computer  code  techniques,  finite 
element  or  finite  difference,  which  divide  both  the  target  and  the  penetrator 
into  a  large  number  of  interconnecting  elements  and  in  each  element  the 
equations  of  motion  are  solved,  along  with  the  equations  of  compatibility  and 
material  constitutive  relations,  as  the  problem  steps  forward  in  time.  A 
large  number  of  codes  are  available  for  both  localized  impact  response  and 
structural  response,  and  there  have  been  several  reviews  of  the  technology 
Iio,  n-131.  The  computer  codes  find  their  niche  in  combination  with 
experiments  in  the  detailed  design  of  structures,  and  as  an  aid  to 
understanding  the  evolution  of  deformation  features  in  impact  events.  By 
themselves  codes  do  not  guarantee  a  complete  numerical  solution  to  a  specific 
impact  event. 

There  has  been  a  keen  interest  in  recent  times  in  one  particular 
class  of  problem  which  is  characterized  by  low  velocity  Impacts  (below  normal 
ordnance  velocities)  of  the  order  of  a  few  hundred  metres  per  second  or  less, 
involving  blunt  or  flat  ended  projectiles  and  relatively  thin  targets  in 
comparison  to  the  projectile  diameter.  such  a  situation  generally  results  in 
the  formation  and  ejection  of  a  plug  from  the  target  and  a  structural  response 
evidenced  by  a  large  permanent  dishing  displacement  of  the  plate.  The 
problem  is  one  of  interest  for  safety  in  many  civil  and  military 
applications.  in  practice  it  is  common  for  the  conditions  to  be  ill-defined 
in  that  only  broad  estimates  of  likely  impact  situations  are  possible,  thus 
the  methods  for  solution  must  be  adaptable  and  efficient.  A  number  of  recent 
papers  (14-23)  present  experimental  data  and  model  solutions  for  cases 
involving  plugging  of  plates  and  both  empirical  approaches  (23)  and  problems 
of  scaling  (241  are  considered. 

in  this  report  two  methods  of  approaching  low  velocity  perforation 
problems  are  presented.  The  first  Involves  assessment  using  a  simple  set  of 
analytic  equations.  The  second  approach  uses  a  structural  model  which 
describes  the  sequence  of  dishing,  stretching  and  plugging  of  the  target  and 
requires  a  computer  solution.  Comparison  with  empirical  data  and  an 
examination  of  the  deformations  occurring  in  impact  events  enables  a 
constructive  approach  to  considerations  of  how  to  handle  such  problems.  It 
is  a  prime  objective  of  this  work  to  provide  a  set  of  simple  tools  which  can 
easily  be  applied,  as  many  of  the  techniques  presented  in  the  literature  are 
in  outline  only  and  require  a  computer  program  to  be  written  by  the  user,  or 
they  require  some  experimentation  which  is  not  always  possible  or  desirable. 


2.  SIMPLE  MODEL  SOLPTIOWS 


If  the  discussion  is  limited  to  ductile  metallic  targets  and 
projectiles  which  are  substantially  non-deforming  during  the  interaction,  then 
it  is  possible  to  derive  simple  equations  for  perforation  depending  on  the 
deformation  mode.  For  a  pointed  projectile  material  is  generally  pushed 
radially  to  the  side  when  the  projectile  impacts  the  target  as  in  Fig. 
l(a).  This  mode  of  failure  is  generally  called  ductile  hole  formation. 

Taylor  til  developed  an  equation  for  the  work  done  in  expanding  a  hole  froei 
zero  diameter,  out  to  the  projectile  diameter  and  this  has  been  shown  (3)  to 
give  a  reasonable  estimate  of  the  work  done  in  penetrating  a  target  by  the 
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ductile  hole  formation  mode.  The  work  done  (W)  la  given  In  terms  of 

projectile  diameter  ( r>) ,  target  thickness  (h)  and  target  strength  (o  )  by 

o 

W  -  f  Djh  o  (l) 

2  O 

If  the  target  Is  thin  (thickness  less  than  the  projectile  diameter) 
then  there  is  a  tendency  for  It  to  bend  forward  m  the  direction  of 
penetration  and  produce  a  dishing  failure.  Pig.  l(b).  In  this  case  work  Is 
done  In  stretching  elements  of  target  out  to  the  diameter  of  the  projectile 
and  also  In  bending  the  plate.  An  equation  for  the  work  done  Is 


W  -  t  Dh  <D  +  «h)o  (2) 

8  O 

For  a  flat  ended  penetrator,  target  material  Is  constrained  to  move 
ahead  of  the  penetrator  so  that  a  plug  is  pushed  out,  Pig.  1(c).  Under  these 
circumstances  the  acceleration  of  material  ahead  of  the  projectile  must  be 
considered  as  well  as  the  plug  shearing  stresses,  and  generally  significant 
thinning  of  the  plug  is  observed.  This  process  has  been  modelled  in  detail 
(6,7, l€— 19 1 ;  however  for  thick  plates  most  of  the  work  done  is  Involved  with 
shearing  of  the  plug.  Thus  an  underestimate  can  be  made  as  the  work  done  in 
simply  shearing  a  plug  of  diameter  D  from  the  plate  of  thickness  h,  namely 


w  -  -1-  Uh2<>  (3) 

For  thin  plates  impacted  by  flat  ended  projectiles  dishing  of  the  plate  will 
also  occur  and  this  is  not  accounted  for  In  equation  (3>. 


The  choice  of  strength  or  flow  stress  parameter  (o  )  to  be  used  In 
equations  (l)  to  (3)  is  determined  by  the  fact  that  the  processes  are  ones  of 
large  strain  deformation  and  dynamic  loading  so  that  a  high  estimate  of 
strength  of  the  material  is  generally  better.  As  most  streas/straln  data  Is 
determined  at  comparatively  low  strain  rates  the  preferred  choices  have  been 

(a)  to  fit  stress/straln  data  to  a  curve  of  the  form 

o  -  o  «R  ( 4 ) 

o 

where  o  and  <  are  stress  and  strain  respectively,  and  the 
constants  o  and  n  are  obtained  using  a  log/log  plot.  The  constant 
oQ  which  is  the  stress  at  a  strain  of  l.o  is  used  for  the  strength. 

(b)  to  use  the  ultimate  tensile  strength  (UTS)  of  the  material 
if  this  is  all  that  is  available. 

using  the  above  methodology  it  is  generally  found  that  equation  (l) 
overestimates  the  work  done  in  penetration,  equation  (3)  vastly  underestimates 
It  for  thin  targets  where  dishing  is  a  major  consideration,  and  equation  (2)  • 

Is  in  the  correct  range.  Examples  of  such  results  are  presented  later  for 

I 
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comparison  with  solutions  using  the  structural  model  below.  Equations  (1)  to 
(3)  have  the  advantage  that  being  simple,  the  solution  can  be  obtained  rapidly 
by  hand. 


3.  STRUCTURAL  MODEL 


Because  the  problem  of  blunt  projectile  impact  on  a  thin  plate 
involves  localized  damage  as  well  as  structural  deformation,  a  model  was 
developed  which  takes  both  aspects  into  account,  and  is  based  on  a  rigid 
plastic  solution  for  the  central  impact  on  an  infinite  beam  developed  by 
symonds  [251.  The  basic  geometry  is  shown  in  Fig.  2.  Impact  of  a 
projectile  of  mass  G  at  a  sufficiently  high  velocity  VQ  onto  an  infinite  beam 
causes  shear  sliding  deformation  at  the  impact  site  if  the  shear  stress  Qp  in 
the  beam  equals  the  shear  yield  stress  of  the  beam  material.  The  beam  bends 
around  a  plastic  hinge  a  distance  z  from  the  impact  site.  Solution  of  the 
equations  of  conservation  of  linear  and  angular  momentum  allows  the  position 
of  the  hinge  and  the  velocity  of  rotation  to  be  obtained  as  a  function  of 
time.  This  geometric  configuration  can  be  transformed  to  that  of  a  plate  if 
appropriate  alterations  are  made  for  the  moment  of  inertia  and  mass 
distributions.  Figure  3  shows  the  geometry  where  a  projectile  of  velocity 
v0,  diameter  2R0,  and  mass  G  Impacts  a  plate  of  thickness  hQ.  The  equations, 
which  have  been  derived  elsewhere  (26!,  are  outlined  below  with  the  method  of 
solution. 


Figure  3  divides  the  process  Into  two  stages.  In  stage  I  the  plug 
and  projectile  slide  at  a  velocity  ▼' ,  greater  than  the  velocity  of  adjacent 
plate  movement,  v.  The  plate  rotates  about  a  hinge  at  position  z.  Equating 
the  impulse  of  the  force,  Qp,  to  the  change  in  momentum  on  either  side  of  the 
plug/plate  interface  gives 

m  R  2 

gv  -  GV’  -  0  0 ■  V’ 

o  2 


m0  -  2»pho 

is  projectile  mass, 
is  projectile  impact  velocity, 
is  the  plug  velocity  (projectile  and  plug  assumed  to  move  together) 
is  the  velocity  of  the  plate  adjacent  to  the  plug, 
is  the  plate  initial  thickness, 
is  the  projectile  radius, 

is  the  hinge  position,  measured  from  the  plug/target  interface, 
is  the  plate  density,  and 
is  time. 


where 


G 

vo 

V' 
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-  f  20  dt  (5a) 

J  P 

-  -f-  (r  +  —  z)v  (5b) 

2  0  3 

(6) 
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equating  the  change  in  angular  momentum  to  the  impulse  of  the  torque 


GV'Z  +  -  m  R  2V'Z  +  7  m  z2V  (R  +  7  z)  -  -2f(M  +  M  )dt  (l> 
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and  MpZ  are  the  plastic  moments  at  the  projectile  radius  and  the 
hinge  position  respectively. 

The  equations  differ  from  those  of  Symonds  [251  in  including  the 
mass  of  the  plug,  in  using  a  mass  distribution  appropriate  to  a  plate  and 
allowing  the  impact  to  operate  over  the  full  projectile  diameter  rather  than 
at  a  central  point. 

The  shear  force  at  the  plug/target  Interface  is  given  by 


Similarly 

gives 


-  gv  z  + 
o 


where 


Q 

-  —  R  h  Y 

(8) 

P 

w| 

o 

and  the  moments  at  the  plug/target  interface  and 
given  respectively  by 

at  the  hinge  position  z  are 

M 

-  7  R  h2Y 

(9a) 

po 

4  0 

H 

-  7  <R  +  z)h  2Y 

(9b) 

pz 

4  0  0 

where  Y  is  the  material  yield  stress  and 

h  is  the  reduced  section  contact  area  at  the  plug/target 
Interface  due  to  shearing. 

Equations  (s)  to  (9)  can  be  combined  to  give  a  quadratic  equation 
for  z  which  can  be  solved  at  each  time  step.  The  initial  position  of  the 
hinge  z  is  found  to  be  R0/3  and  this  increases  with  time.  solution  of  the 
equations  shows  that  the  velocity  of  the  plate,  V,  Increases  with  time  until 
it  either  shears  out  as  a  plug  or  it  equals  the  plug/projectile  velocity  V' . 

At  this  point  shearing  of  the  plug  ceases  as  the  relative  plug/plate  velocity 
is  zero  and  the  plate  continues  to  stretch  as  a  membrane  with  bending 
continuing  to  occur  at  the  hinge  position  z.  This  is  referred  to  as  stage  II 
and  the  conditions  are  shown  in  Fig.  3. 


To  simplify  the  calculations  in  stage  II  it  was  assumed  that  at  the 
projectile  radius  the  yield  stress  in  tension  is  exceeded  across  the  full 


section  of  the  plate. 


The  membrane  force  Fp  is  then 


F  «  «R  h  Y  <10) 

p  O 
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and  the  bending  moment  at  z  -  o  Is  zero.  At  the  hinge  position,  z,  the  full 
plastic  moment,  MpZ  is  again  given  by  equation  (9b).  However,  because  of  the 
membrane  force  F„  a  reduced  moment  M'  acts  (27,261,  which  is  given  in  terms 
of  V  as  P  pz 


pz 


R  h 

M  (2  -  (T— 2 - r— )2> 

pz  (r  +  z)h 


(li) 


The  equations  for  angular  and  linear  momentum,  respectively,  in 
stage  II  are 


■7  Gzav  Cos  •  +  4  m  R  2zAV  Cos  e  +  \  m  z2  (R  +  7  z)  IV  Cos  # 
2  400  £004 


-  -  (M-  +  F  z  sin  e)  At 

P  P 

GAV  +  77  m  R  ZAV  +  7  m  z  (R  +  7  z)  AV  +  a  (r  +  7  z)  vaz  »  0 
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(12) 

(13) 


where  e  is  the  angle  through  which  the  plate  is  bent. 
At  is  the  time  step, 
av  is  the  change  in  velocity  and 
az  the  change  in  hinge  position. 


Terms  of  the  type  AZAV  have  been  ignored.  As  z  is  known  at  the 
start  of  stage  II,  equation  (12)  can  be  solved  for  AV  and  equation  (13) 
for  az  at  each  time  step.  Thus  at  each  time  step  the  moments,  change  in 
velocity,  change  in  hinge  position  and  the  new  angle  «  are  calculated.  As 

the  plate  is  deforming  as  a  membrane  in  tension  and  as  there  is  a  reduced 
section  at  the  plug/plate  interface,  the  expected  failure  is  a  tensile 
fracture  at  that  point.  A  suitable  failure  criterion  is  required  to 
calculate  the  amount  of  membrane  stretching  before  f  acture. 


Consider  the  simplification  that  in  stretching  as  a  membrane  in 
tension  the  plate  behaves  as  a  linearly  tapered  tensile  sample  and  has  linear 
strain  hardening  of  the  form 


o  -  «  +  ft 


(14) 


where  o  is  the  flow  stress  and 
<  is  engineering  strain, 
a  and  t  are  constants. 

The  situation  is  depicted  in  Fig.  4,  and  considering  the  plate  as  stretching 
radially  from  the  plug  diameter,  the  angle  of  taper  is  2*.  The  strain  at 
any  position,  under  these  conditions,  is  inversely  proportional  to  radial 
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position,  and  the  mean  strain,  c  ,  in  the  membrane  can  be  related  to  the 
strain  at  the  inside  radius  RQ  and  the  maximum  radius  of  yielding,  Rj ,  by 
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R,  R 
1  o 


(R,  -  R  r 
1  o 


+  n 


(IS) 


In  the  particular  case  of  linear  work  hardening 


R 
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Thus  the  mean  strain  is  given  in  terms  of  the  strain  at  RQ  by 


a  (  a  +  g 


a  +  g 1 


-  (- 


in  (- 


1>-V 


(17) 


At  fracture  equals  (f,  the  fracture  strain  in  tension  of  the 
material.  Thus  equation  (17)  gives  a  value  for  the  mean  membrane  strain  in 
the  dished  material  at  fracture  in  terms  of  the  material  characteristics 
a,  g  and  if.  «f  can  be  obtained  from  a  simple  tension  test  and  a  and  g  by 
fitting  the  material  stress/straln  characteristics  to  equation  (is). 

In  carrying  out  the  calculations  of  stage  II,  the  change  in  angle 
of  dishing  as  in  each  step  is  converted  into  an  increment  of  mean  membrane 
strain  a<  .  At  each  step  the  emulative  mean  membrane  strain  within  the 
dished  region,  which  extends  to  the  hinge  at  position  z,  is  summed  till 
equation  (17)  is  satisfied  with  «  equal  to  »  for  the  material.  At  this 
point  the  plate  is  considered  per?orated  due  to  fracture  at  the  periphery  of 
the  plug.  This  approach  to  failure  is  very  simplified  in  many  aspects, 
including  the  use  of  engineering  strains,  the  assumption  of  a  uniaxial  stress 
state  in  a  tapered  tensile  geometry  for  the  membrane  stretching  of  the  plate, 
the  incrementing  of  strains  with  changes  in  dishing  angle,  and  the  simplified 
approach  to  the  material  properties.  The  method  has  however  provided  a 
consistent  picture  when  applied  to  a  number  of  cases  which  have  been  examined 
and  is  used  at  this  stage  for  want  of  a  more  rigorous  incremental  procedure. 


As  the  model  is  for  the  consideration  of  impacts  below  conventional 
ordnance  velocities  consideration  must  be  given  to  the  possible  effects  of 
clamping  of  the  target  plate  if  this  is  close  enough  to  the  impact  site  to 
Influence  the  plate  deformation.  This  is  easily  done  computationally  by 
including  an  option  for  the  maximum  value  of  the  hinge  position  z  as  the 
position  of  the  clamp. 

The  program  listing  of  the  structural  model  called  DASB  is  given  in 
the  Appendix  with  a  typical  listing  of  input  data.  A  variety  of  output 
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parameters  are  provided  for,  and  a  typical  output  listing  is  given  for  the 
input  data  provided. 


«.  DISCUSSION  AMP  COMPARISON  OF  MODEL  RESULTS 


The  simple  analytical  models  presented  restrict  the  deformation 
considerations  to  one  simplified  mode  in  each  case.  The  structural  model  on 
the  other  hand  allows  for  bending  and  stretching  of  the  plate  as  well  as 
shearing  of  a  plug.  The  mechanism  of  plug  separation  for  thin  plates 
described  in  the  model  can  be  by  a  tensile  failure  or  a  shear  fracture 
depending  on  the  material  and  impact  velocity,  and  this  is  in  accord  with 
observations  of  plug  separation  in  thin  plates.  For  example.  Fig.  5  taken 
from  Levy  and  Goldsmith  1211  is  one  of  many  examples  clearly  illustrating  the 
tensile  separation  of  a  plug.  Shear  failures  by  sliding  off  the  plug  are 
described  in  the  model  for  thick  targets  and  higher  velocities  of  Impact  and 
this  occurs  directly  from  stage  X.  Thus  the  structural  model  conforms 
qualitatively  with  experience. 

Shadbolt  et  al.  [l«l  present  a  comparison  of  several  models  with 
experimental  data  for  perforation  velocity  as  a  function  of  thickness  for  mild 
steel,  aluminium  alloy  and  stainless  steel  targets.  The  same  data  are 
compared  with  the  analytical  solutions  for  ductile  hole  formation  failure, 
equation  (l),  dishing  failure,  equation  (2>,  and  plugging  failure,  equation 
(3),  as  well  as  the  structural  model  in  Fig.  6.  It  is  noted  that  in  each 
case  the  simple  plugging  solution  greatly  underestimates  the  critical 
velocity,  whereas  the  experimental  data  generally  lie  between  the  ductile  hole 
formation  and  dishing  solutions.  The  structural  model  gives  similar  results 
to  the  dishing  solution.  Examination  of  the  energy  absorbed  by  different 
deformation  modes,  membrane,  shearing  and  bending  work,  in  the  structural 
model  gave  very  similar  results  to  the  approximate  experimental  estimates  of 
Corran  et  al.  [151  for  the  mild  steel  shots.  Whilst  the  use  of  the 
structural  model  is  a  better  approach  these  comparisons  indicate  the  degree  to 
which  the  simpler  solutions  can  be  used  as  useful  approximations. 


Scale  up  to  larger  size  missiles  can  be  assessed  to  some  degree 
using  the  results  of  Neilson  123)  produced  at  the  Atomic  Energy  Establishment, 
Winfrith,  UK,  for  impacts  on  mild  steel  plates  from  1  to  25  mm  thick. 

Missile  diameters  up  to  85  mm  and  masses  to  20  kg  were  used.  Table  I 
compares  computations  using  the  structural  model  with  the  experimental 
data.  Given  the  need  to  estimate  material  properties  the  agreement  is 
good.  For  the  tests  A  and  B  using  the  lower  strength  steel  the  only 
difference  was  the  span  width  of  the  target  which  Influenced  the  model 
solution  in  the  correct  manner.  For  all  other  examples  the  span  was 
sufficient  that  it  did  not  Influence  the  calculated  critical  perforation 
velocity. 


As  the  objective  is  to  give  a  more  general  solution  which  can  be 
applied  to  estimate  behavior  with  limited  knowledge  of  material  data  and  only 
a  broad  appreciation  of  the  likely  impact  conditions  the  structural  model  was 
also  tested  for  the  impact  of  spheres  on  hollow  steel  tubes.  Photographs  of 
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typical  failures  observed  by  Xlaoqlng  and  stronge  (22)  showed  that  the  basic 
deformation  inodes  involve  dishing  and  plugging,  so  that  even  though  there  is 
curvature  in  the  target  it  is  worth  examining  the  structural  model  to  see  if 
it  can  be  used  to  estimate  ballistic  limit  velocities  for  various  tube  wall 
thicknesses  and  ball  diameters.  The  comparison  of  experimental  data  with  the 
calculated  solutions  in  Fig.  7  again  shows  a  correct  order  of  magnitude 
estimate  and  correct  trends  with  projectile  diameter  and  tube  thickness. 


The  diverse  cases  examined  show  the  areas  of  applicability  of  the 
structural  model  and  that  it  can  be  relied  upon  to  give  useful  broad 
predictions  of  behaviour  within  what  is  acceptably  possible  from  the  theory  of 
plasticity.  The  structural  model  and  the  simple  analytical  solutions  are  in 
a  form  which  can  easily  be  applied  by  practising  engineers.  A  particularly 
valuable  aspect  is  that  rapid  estimates  of  the  influence  on  design  variations 
of  factors  such  as  projectile  diameter  and  mass  and  target  thickness  and 
material  characteristics,  can  be  obtained. 


5 .  CONCLUSION 


This  report  outlines  two  approaches  to  making  estimates  of  the 
resistance  of  thin  plates  to  perforation  by  low  velocity  impacts.  one 
technique  uses  a  set  of  analytical  equations  and  the  alternative  method  uses  a 
structural  model  which  requires  a  computer  solution.  The  approaches  are 
compared  with  published  experimental  data  for  projectile  impacts  to  indicate 
the  level  of  certainty  expected  from  predictions.  The  essential  purpose  of 
the  methods  is  to  give  guidance  on  the  overall  level  of  protection  offered  by 
a  plate  where  the  material  properties  are  not  accurately  known,  impact 
conditions  are  only  broadly  defined  and  experiments  are  not  possible. 
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TABLE  I 


Comparison  o£  Neilson  (231  Data  with  Structural  Model  Calculations 


SPAN 

DESIGNATION  WIDTH 
(mm) 

TAROT 

THICKNESS 

(ran) 

IRQJECTIIE 

DIAMETER 

(ran) 

PROJECTIZE 

MASS 

(kg) 

OTS* 

(MPa) 

vc 

EXPEROENT 

(ms-1) 

VC 

CAL3HAXICV 

(ms-1) 

A 

220 

1.0 

32 

4.3 

320 

12.5 

17.5 

B 

127 

1.0 

32 

4.3 

320 

9.5 

12.5 

C 

S8S 

15 

63 

19.8 

470 

80 

75 

D 

588 

25 

63 

19.8 

470 

102 

101 

E 

1500 

6 

85 

9.72 

473 

79 

65 

F 

1500 

12 

85 

9.72 

473 

99 

91 

G 

1500 

6 

43 

1.22 

473 

98 

91 

H 

280 

1.0 

40 

1.0 

420 

43 

42 

X 

280 

3.0 

40 

1.0 

420 

83 

67 

*  Parameters  used  for  structural  model  calculations  (equation  14) 


A  and  B 

Yield 

strength 

<Y) 

-  o  -  220 

MPa, 

P 

-  360 

MPa,  <f  - 

.37 

C  to  I 

Yield 

strength 

(Y> 

■  a  ■  340 

MPa, 

P 

>  260 

MPa,  t f  - 

.37 

APPEKDIT 


The  computer  program  for  solution  of  the  structural  model  Is  called 
DAS a .  The  Input  file  Is  called  DASI.DAT  and  there  are  three  output  files 

DASP.DAT,  DASA.DAT  and  DASO.DAT.  The  units  for  Input  data  are  Indicated  In 
Table  Ai  and  all  calculations  and  output  are  In  SI  units.  There  are  a 
limited  number  of  comment  statements  In  the  program  to  Identify  the  major 
steps  and  the  output  parameters.  For  problems  where  the  width  of  span  Is 
unimportant,  99999.  will  assume  the  plate  Is  of  Infinite  span.  In  other 
cases  the  maximum  hinge  position  Is  limited  by  the  Input  width  of  span.  Table 
A2  gives  a  typical  output  listing  for  the  data  of  Table  Al. 


TABLE  Al 

Typical  Input  -  DASI.DAT  (Mild  steel  target) 


Parameter 

▼alue 

Units 

Projectile  diameter 

12.7 

mm 

Target  thickness 

3.0 

mn 

Target  density 

7.8 

g  cm-3 

Target  yield  strength 

220. 

MPa 

Projectile  mass 

.  0  346 

kg 

Impact  velocity 

100. 

ms-1 

o  (equation  14) 

220. 

MPa 

t  ( equation  1 4 ) 

360. 

MPa 

Fracture  strain  ( <  f ) 

.37 

dimensionless 

span  width 

99999. 

mm 

TABLE  A2 

Typical  Output  -  DASO.DAT* 


T  H  Z  V  X 

0.000118  0.00202  0.0184  35. C  0.00C7 


PEE 

WS 

WB 

WT 

REE 

TE 

23.9 

12.5 

62.6 

30.3 

7.1 

136.4 

*  see  program  listing  for  identification  of  symbols 


o  o  o  o  o 


LISTING  OF  PROGRAM  CASH 


CASH. FOR 

PROGRAM  FOR  PENETRATION  OF  THIN  TARGETS 
WRITTEN  BY  RAYMOND  L  WOODWARD, DEPT  OF  DEFENCE, 

MATERIALS  RESEARCH  LABORATORIES, 

MARIBYRNONG  1986 

OPEN ( UNIT-1,  FILE- 'CASI .CAT’ ,STATUS-'OLD' ) 

OPEN(  UNIT-2, FILE- 'DASF.DAT' , STATUS- 'NEW' ) 

OPEN  (UNIT-3,  FILE- 'CASA. DAT' ,STATUS-'NEW' ) 

OPEN (UNIT-4,  FILE- 'DASO.DAT' ,STATUS-'NEW' ) 

21  FORMAT  ( F10 . 2/F10 . 2/F6 . 3/F7 . 1/FI  0 . 3/T7 . 1/F7 . 1/F7 . 1/F5 . 3/F8 . 1 ) 

23  FORMAT(F10.8,4X,F10.1,4X,F10.1,4X,F10.1,4X,F10.1,4X,F8.5,4X,F12.3 
1,4X,F10.5,4X,F10.5) 

22  FORMAT(F10.8,4X,F10.5,4X,F10.5,4X,F10.5,4X,F10.1,4X,F10.1,4X,F6.1 
1,4X,F6.1) 

24  FORMAT)  3X,F10.6/3X,F10.6/3X,F8.1/3X,F12.0/F10.5/3X,F6.1, 

1/3X,  F12 . 0/3X,  F12 . 0/3X,  F6 . 3/3X,  F8 . 1/3X,  F12 . 1 ) 

25  FORMAT(2X,F10.6,2X,F8.5,2X,F8.4,2X,F6.1,2X,2X,F8.4,F10.1,2X,F10.1 
1,2X,F10.1,2X,F10.1,2X,F10.1,2X,F10.1,2X,F10.1,2X,F10.1) 

28  FORMAT(F10.8,4X,F10.5,4X,F10.5,4X,F6.2,4X,F10.1,4X,F10.1,4X,F10.1 
1 , 4X, F10 . 1 , 4X, F10 . 1 , 4X, F6 . 1 ) 

26  FORMAT (T4, 'PROJECTILE  THROUGH') 

27  FORMAT) T4 , ' PROJECTILE  STOPPED') 

29  FORMAT) T4, 'TIME  EXCEEDED' ) 

64  FORMAT(6X, 'T'lOX, 'H'lOX, 'Z'lOX, 'V'lOX, 'X'9X, 'PKE'9X, 'WS'9X, ' 
1WB'9X'WT'9X,  'RKE'9X,  'TE' ) 

65  FORMAT) 6X, 'T'12X, 'Z'12X, 'X'14X, 'XP'12X, 'PKE'lOX, 'TE'lOX,  'V' 

110X, 'VP' ) 

66  FORMAT) 6X, 'T'14X, 'PKE'12X,  'WS'llX, 'WB'12X, 'RKE'9X, 'H' 

114X, 'BMI'llX, 'Zl' ,11X, 'Z2' ) 

67  FORMAT) 6X,  'T'12X,  'Z'12X,  'X'12X,  'TH0'12X,  'WB'12X,  'RKE'12X,  'WT' 
111X,'PKE'12X,'TE') 

READ( 1,21 )DO,HO,RHO,SO,ASS,VEL,ALPA,BETA,EF,SPAN 
C 

C  DO  PROJ  DIAM  MM 

C  HO  TARGETT  THICKNESS  MM 

C  RHO  TARGET  DENSITY  GM/CC 

C  SO  TARGET  YIELD  STRENGTH  MPA 

C  ASS  PROJ  MASS  KG 

C  VEL  PROJ  VEL  M/S 

C  ALPA  STRAIN  HARDENING  RELATION  TARGET 

C  BETA  YS-ALPA+BETA*STRAIN  (MPA) 

C  EF  FRACTURE  STRAIN  IN  TENSION  TEST 

C  SPAN  TOTAL  WIDTH  OF  TEST  PIECE  MM  IF  UNIMP.  USE  99999. 

C 

HO-BO/IOOO 

RD-DO/2000 

RHO-RBO*1000 

80-80*1000000. 

ALPAmALPA*1000000 . 

BETA-BETA*1000000 . 

ZMAX-) SPAN-DO  )/2000 

EMF-EF*  ( (ALPA* (ALPA+BETA*EF )/( ( BETA*EF ) **2 ) )  *ALOG( 1+BETA*EF/ 
1ALPA) -ALPA/ ( BETA* EF) ) 

SKE-ASS* ( VEL**2  )/2 
TWO 


noon  no 


DO-6.28318*RHO*HO 
C  QI  IS  SHEAR  FORCE  INTEGRAL 
01-0 

C  K  IS  BOOIN5  MOMENT  INTEGRAL 
BMI-0 
NB-0 
MS-0 
MT-0 

F1-*SS*VEL 

RD2-RO**2 

R03-*0**3 

F2-ASS+O40*  ( RO**2 )/2 

X— 0 

XP-0 

Zl-0 

Z2-0 

C  TIME  STEP  FOR  FIRST  STAGE 

DT-.008019*ASS*VEL*RHO*RO/(  (ASS+4.044*RBD*BO*(RO**2) )*SO) 

WRITES  OUT  INPUT,  SKE- INITIAL  KE 

WRITE( 4,24 )DO,HO,HHO,SO, ASS, VEL,ALPA,BETA,EF, SPAN, SKE 
WRITE(4,64) 

HRITE(2,65) 

WRITE(3,66) 

INITIAL  PHASE  OF  BOOING  AND  SHEARING  TILL  VEL  BOOING  PLATE 
EQUALS  PLUG/PROJ  VEL 

DO  72  1-1,500 
QP-1.814*H*SO*RO 
QI-QI+2*QP*DT 
VP-( Fl-QI )/T2 
IF(VP.LE.V)00  TO  73 
A— 3*QI 
B— RD*QI-2*BMI 
C-6*BMI*RO 

Zl-(-B+SQRT( (B**2)-4*A*C)  )/(2*A) 

Z2-(-B-SQRT( (B**2)-4*A*C)  )/(2*A) 

Z-Z2 

V-2*QI/(  DD*Z*  ( KHZ/3 ) ) 

X-X+V*DT 
XP-XF+VP*DT 
H-BCMXP-X) 

PKB-F2*(VP**2)/2 
WS-WS+2*QP*  (VP-V) *DT 
IF(Z.GT.O)GO  TO  32 
WB-MB+0 
GO  TO  33 

32  MB-MB+l .  5708*SO*V*DT*  ( RO*  ( H**2 )+ ( KHZ )  *  ( H0**2 )  )/l 

33  RKE-1 . 0471*HO*RHO*Z* ( ROfZ/4 ) * ( V**2 ) 

TE-PKE4MS+WB+RKE 
T-T+DT 

BMI-BMI+(RO*  ( H**2  )♦(  KHZ )  *  ( H0**2 ) )  * .  7853*SO*DT 
C 

c  catrvr  msf.mt 

C  T-TIME , Z-HINGE  POSITION, X-DISPL.  PLATE , XP-DISPL .  PLUG, 

C  PKE-PROJ.  KE,TE- TOTAL  e«GT,V-VEL  PLATE, VP-VEL  PIUG. 

WRITE(2,22)T,Z,X,XP,PKE,TE,V,VP 

C 


ononn  non  oooo 


C  OUTPUT  CASA.  CAT 

C  T-TIME , PKE-PROJ .  KE,WS-WORK  SHEAR, WB-WORK  BENDING , RKE-ROTATIGNAL  KE, 

C  H- PLUG/PLATE  CCNTACT  LENGTH,  BMI-BB4D  MTT  INTEGRAL, 

C  Zl-INVALID  QUADRATIC  FOOT.Z2-VALID  BOOT  HINGE  POSITION. 

NRITE(3,23)T,PKE,WS,NB,RKE,H,BMI,Z1,Z2 
ir(H.LE.O)GO  TO  75 
IF(VP.LT..l)GO  TO  76 

72  CCNTINUE 

SECOND  STAGE  OF  BENDING  AND  STRETCHING  TILL  FRACTURE 
CONDITION  EXCEEDED 

73  V-VP 
THR-ATAN2(X,Z) 

ENTR-Z/OOS(THR) 

FT-6.2832*RO*H*SO 

THB-THR 

THP-THR 

EREF-( 1/COS (THR) )-l 

TIME  STEP  SET  FOR  SECOND  PHASE 

DTWHO/(25*VEL) 

EM-0 

WRITE (2, 67) 

D052  EM.,  5000 

SNP-1 . 5708*SO* ( ROfZ ) * ( HO*  *2 ) 

SMPP-SMP* ( 1- ( RO*H/ ( (RO+Z)*HO) )**2) 

DV-( -SMPP-FT*Z*SIN(THE)  )*DT*2/<  <F2+DK)*Z*(R04-Z/4)/3)*Z*C0S(THE) ) 

DZ-(  -F2*DV-EMO*DV*Z*  ( RO+Z/3 )  *DV/2 )/( 6 . 283*  ( RD* .  6667*Z )  *HO*RHO*V) 

V-V+DV 

X-X+V*DT 

THE-AXAN2(X,Z) 

DEN- 1/COS  (THE) -l/COS  (THP) 

WB-WB+( SKPP+FT*Z*SIN(THE ) )*(THE-THP) 

Z-Z+DZ 

IF( SPAN. GT. 99998. )GO  TO  109 
IF(Z.LE.ZMAX)00  TO  109 
Z-ZMAX 

L09  PKB-F2*(V**2)/2 

RKE-1 . 0471*HO*RHO*Z*  ( HOfZ/4 )  *  ( V**2 )/( COS(THE)  **2 ) 

WTWNT4-FT*O0S  ( THE )  *  ( ZP/COS  ( THE ) -ZP/COS  ( THP ) ) 

ZP-Z 

THFtAXAN2(X,Z) 

TE-PKE+WB*WB+OXE4WT 

T-T+DT 

EM-DHDOI 

TBO-THE*57.296 

IF(EHF.LE. (EM-EREF) )GO  TO  75 

IF(V.LT. .1)00  TO  76 

DO  83  1-1,200 

INCVI*25 

IF(N- INC >52,53,83 
83  CONTINUE 

OUTPUT  DASF.DAT 

T-TIME, Z-HINGE  POS,X-DISPL.  PLATE, THO-BEND  ANGLE  PLATE  DEGREES 
WB-WORX  BENDING, RKE-ROTAX.  KE,WT-ME»RANE  WORK, 


nnnnnnn 


C  PKE-PROJ.  FE.TE- TOTAL  ENERGY. 

S3  WRITE(2,28)T,Z,X,TH0,WB,RKE.WT,PKE,TE 
52  CONTINUE 
WRITE(4,29) 

00  TO  77 

75  WRITE(4,26) 

00  TO  77 

76  WRITE(4,27) 

OUTPUT  DASO.QAT 

MAJOR  OUTPUT  PARAMETERS, T-TIME,  H-PLUG/EARGET  CONTACT  LENGTH, 
Z-FINAL  HINGE  POSITION,  V-PLUG/PROJ  VEL,  X-DISHING  AMT, 
PKE-PROJ  RE,  WS-WCRK  SHEAR,  WB-MORK  BENDING,  WT-TENSION  WORK, 
RKE- ROTARY  RE  OF  PLATE,  TE-TOTAL  ENERGY 


77  WRITE(4,25)T,H,Z,V,X,PKE,WS,WB,WT,RKE,TE 
STOP 
END 


FIGURE  1 


Mechanisms  of  failure  of  ductile  metal  targets. 

(a)  Ductile  hole  formation  failure,  (b)  dishing  failure  and 
(c)  plugging  failure. 


FIGURE  2 


Impact  of  a  mass  G  with  velocity  v0  on  an  Infinite  beast. 
Shearing  and  bending  occurs  at  the  impact  sits  and  bonding 
at  a  hinge,  distance  s  from  Impact. 


FIGURE  3 


FIGURE  4 


t 


Impact  of  a  projectile  on  an  Infinite  plate. 

(a)  schematic  section,  (b)  stage  I  in  which  shearing  of  the 
plug  occurs,  and  (c)  Stage  II  in  which  membrane  stretching 
of  the  plate  occurs. 


Treatment  of  the  plate  In  stage  n  as  a  tapered  tensile 
sample  with  material  properties  characterised  as 
rigld/llnear  work  hardening.  The  arrows  indicate  the 
application  of  the  tensile  force  and,  as  the  radial 
stretching  of  a  plate  le  being  considered,  the  angle  of 
taper  is  2*. 


I  «  •  I 

t«cumU  >»m> 
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FIGURE  <  comparison  of  Solution*  using  th#  modal*  for  ductllo  hoi* 

formation  (DBF) ,  aquation  (1),  dishing,  aquation  (2), 
plugging,  aquation  (3>,  and  th*  structural  modal  with  th* 
axparlmantal  data  of  Corran  at  al.  (IS, 1(1.  (a)  Mild  ataal 

targets,  (b)  aluminium  alloy  targat*  and  (c)  atalnia** 
steal  targets. 
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figure  7  Comparison  of  empirical  data  from  Xlaoqing  and  Strongs  122] 

for  tubes  of  different  wall  thickness  and  for  different 
projectile  diameters  with  predictions  using  the  structural 
model.  Results  are  offset  at  each  missile  diameter  to 
distinguish  the  different  tube  size  results  more  clearly. 


